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Abstract. We begin by reviewing the theory of thermal, neutrino-driven proto-neutron 
star (PNS) winds. Including the effects of magnetic fields and rotation, we then derive 
the mass and energy loss from magnetically-driven PNS winds for both relativistic and 
non-relativistic outflows, including important multi-dimensional considerations. With these 
simple analytic scalings we argue that proto-magnetars born with ~ millisecond rotation 
periods produce relativistic winds just a few seconds after core collapse with luminosi- 
ties, timescales, mass-loading, and internal shock efficiencies favorable for producing long- 
duration gamma-ray bursts. 
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1. NEUTRINO-DRIVEN PNS WINDS 

After a successful core-collapse supernova (SN), a hot proto-neutron star (PNS) 
cools and deleptonizes, releasing the majority of its gravitational binding energy 
(~ 3 x 10 53 ergs) in neutrinos. With initial core temperature T > 10 MeV, a PNS 
is born optically-thick to neutrinos of all flavors because the relevant neutrino- 
matter cross sections scale as a un oc e 2 oc T 2 , where e v is a typical neutrino energy. 
Indeed, because neutrinos are trapped, a PNS's neutrino luminosity L v remains 
substantial and quasi-thermal for a time after bounce tkh ~ 10 — 100 s, as roughly 
verified by the 19 neutrinos detected from SN1987A 20 years ago [1],[2]. Although 
this Kelvin-Helmholtz (KH) cooling epoch is short compared to the time required 
for the shock, once successful and moving outward at ~ 10 4 km/s, to traverse 
the progenitor stellar mantle, tkh is still significantly longer than the time over 
which the initial explosion must be successful. While the specific shock launching 
mechanism is presently unknown, it must occur in a time i< 1 s <C tkh after 
bounce for the PNS to avoid accreting too much matter. 

Thus, even after the SN shock has cleared a cavity of relatively low density mate- 
rial around the PNS, L u remains substantial. Detailed PNS cooling calculations [3] 
show that the electron neutrino (antineutrino) luminosity L Ue ( De ) is ~ 10 52 erg/s at 
t ~ 1 s and declines as oc t~ x until t ~ t K h, after which L v c D \ decreases exponentially 
as the PNS becomes optically thin. This persistent neutrino flux F v jf> e \ continues 
to heat the PNS atmosphere, primarily through electron neutrino (antineutrino) 



absorption on nuclei iv e + n — > p + e~ and v e +p — > n + e + ). Because the inverse, 
pair capture rates dominate the cooling, which declines rapidly with temperature 
(q~ oc T 6 ) and hence with spherical radius r, a region of significant net positive 
heating (q = q + — q~ > 0) develops above the neutrinosphere radius R v . This heat- 
ing drives mass- loss from the PNS in the form of a thermally-driven wind [4]. To 
estimate the dependence of the resultant mass-loss rate (M t h) on the PNS proper- 
ties explicitly, consider that in steady state the change in gravitational potential 
required for a unit mass element to escape the PNS (GM/R U ) must be provided 
by the total heating it receives accelerating outwards from the PNS surface: 

GM dr 
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where M is the PNS mass, v r is the outward wind velocity, and q is per unit mass. 
Because q is quickly dominated by heating from neutrino absorption, which scales 
as q + oc F Ue a nu oc L^e^/Anr 2 , we see that equation (CQ) implies that 

GM L^el too L Ue el e 
ti v M t h JRu M t h 

where we have used M t ^ = 4irpr 2 v r for a spherical wind, p is the mass density, H 
is the PNS's density scale height, crudely defines a mean electron neutrino or 
antineutrino energy, and a subscript V" denotes evaluation near R u . Neglecting 
rotational support and assuming that the thermal pressure P is dominated by 
photons and relativistic pairs (which also becomes an excellent approximation 
as the density plummets abruptly above the PNS surface), we have that H u ~ 
Pu/PuQu oc T^Rl/Mp u , where g v oc M/R 2 V is the PNS surface gravity and T v oc 
(L Ue el e / Rl) 1 / 6 is the PNS surface temperature. T u is set by the balance between 
heating and cooling at the PNS surface (T^ oc if = q + oc L Vii e 2 /R 2 ). Inserting these 
results into equation (j2J) and including the correct normalization from the relevant 
weak cross sections, one finds the expression for M t h first obtained by ref [4]: 

Mth ~ 10" 4 L^ 3 e^M^Rf M /s, (3) 

where L 52 = L Ve x 10 52 erg/s, e 10 = 10e„ e MeV, R v = 10R 1Q km, and M = 1AM 1A M Q . 

Endowed with an enormous gravitational binding energy and a means, through 
this neutrino-driven outflow, for communicating a fraction of this energy to the 
outgoing shock, a newly-born PNS seems capable of affecting the properties of 
the SN that we observe. However, a purely thermal, neutrino-driven PNS wind 
is only accelerated to an asymptotic speed of order the surface sound speed: 
~ c S)i/ ~ ^2kT u /m p w O.lL^^e^-Rio 1 ^ 6 c. Thus, the efficiency r\ relating wind 
power i? t h ~ M t h(f^) 2 /2 to total neutrino luminosity [L v ~ 6L Ve ) is quite low: 

_ E th 5 5/ 6 n /3 4/3 2 

V = -j 10 U2 ClO #10 M 1.4- ( 4 ) 



In particular, although neutrino energy deposited in a similar manner may be 
responsible for initiating the SN explosion itself at early times (i.e., the neutrino 
SN mechanism [5]), rj drops rapidly as the PNS cools. Quasi-spherical winds of this 
type are therefore not expected to affect the SN's nucleosynthesis or morphology 
(although the wind itself is considered a promising r-process source [4]). 

2. MAGNETICALLY-DRIVEN PNS WINDS 

Some PNSs may possess a more readily extractable form of energy in rotation. 
A PNS born with a period P = P ms ms is endowed with a rotational energy 
E rot ~ 2 x 10 52 P~ s 2 i? 2 M 14 ergs, which, for P < 4 ms, exceeds the energy of a typical 
SN shock (~ 10 51 ergs). Given a mass loss rate M and torquing lever arm u T , a wind 
extracts angular momentum J from the PNS at a rate J ~ Quj^M, where Q = 27i/P 
is the PNS rotation rate. With the PNS's radius R v as a lever arm and the modest 
thermally-driven mass- loss rate given by equation ([3]), the timescale for removal of 
the PNS's rotational energy, tj = Jj J ~ MR 2 ,/ Mu 2 ~ Mj Mth, is much longer than 
tkh- However, if the PNS is rapidly rotating and possesses a dynamically-important 
poloidal magnetic field B p (through either flux-freezing or generated via dynamo 
action [6]), then both M and u T can be substantially increased; this reduces tj, 
allowing efficient extraction of E TOt . 

For magnetized winds u T is the Alfven radius u>a, defined as the cylindrical 
radius where pv 2 /2 first exceeds B 2 /8n [7]. The magnetosphere of a PNS is 
most likely dominated by its dipole component, with a total (positive-definite) 
surface magnetic flux given by $b = 2ttB v R 2 , where B v is the polar surface field. 
To estimate uja for magnetized PNS outflows recognize that mass and angular 
momentum are primarily extracted from a PNS along open magnetic flux. For 
an axisymmetric dipole rotator this represents only a fraction pa 2(7r#L CFL )/47r — 
R v /2ujy of $B ; where 6*lcfl ~ \J Rv/^y is the latitude (measured from the pole) 
at the PNS surface of the last closed field line (LCFL), uoy is the radius where the 
LCFL intersects the equator (the "Y point"), and we have assumed that ojy ^> R v 
(#lcfl Cl). Plasma necessarily threads a PNS's closed magnetosphere and cannot 
be forced to corotate superluminally; thus ujy cannot exceed the light cylinder 
radius u L = c/Q = 48P ms km, making it useful to write the PNS magnetosphere 's 
total open magnetic flux as $B,o P en ~ ^B v R 2 v {R v jijJi){ujy /c^l)" 1 • Now, the overall 
latitudinal structure of a PNS magnetosphere (i.e., the allocation of open and closed 
magnetic flux, and the value of wy/wl) is primarily dominated by the dipolar 
closed zone. However, recent numerical simulations [8] show that where the field 
is open it behaves as a "split monopole". In this case the poloidal field scales as 
Bp ~ $B,o P cn/^ 2 ~ 0.2B U P~^ R 10 (lu y / lui,)' 1 (R u /r) 2 , rather than the dipole scaling 
oc (R u /r) 3 . The constant of proportionality is chosen to assure that B P (R U ) B v in 
the limit of vanishing closed zone {u^uiy — > Ry) and is in agreement with numerical 
results (see eq. [28] of ref [8]). 



2.1. Non-Relativist ic Winds and Asymmetric Supernovae 



Non-relativistic (NR) magnetically-driven winds reach an equipartition between 
kinetic and magnetic energy outside oja such that the kinetic energy flux at uja 
(Mv r {u)A) 2 /2) carries a sizeable fraction of the rotational energy loss extracted by 
the wind's surface torque E TOt = Jfi = MQ 2 u\; thus, we have that v t {ua) ~ Quja- 
Combining this with the modified monopole scaling for B p motivated above and 
mass conservation Mq = pr 2 v r (Mq is the mass flux per solid angle) we find that: 

u A /R v ~ B{i*P^M^R%\uy/^)-\ (5) 

where Mq = Mq_± x 1O~ 4 M0S -I sr -1 , B v = P 15 x 10 15 G, and we have concentrated 
on the open magnetic flux that emerges nearest the closed zone (polar latitude 
~ #lcfl) and which thereby dominates the spin-down torque. 

From equation (J5J) we see that winds from rapidly rotating PNSs with 
surface magnetic fields typical of Galactic "magnetars" (B v ~ 10 14 — 10 15 
G) possess enhanced lever arms for extracting rotational energy [9]. Fur- 
thermore, their total outflow power E™ « E rot « 2n8 2 jCFh MQfl 2 u' A « 

10 49 Bis P^^Mq^rII ' 3 (o;y/^l) _3 ergs/s dominates thermal acceleration 

(E™>E th ) for'p 15 > ; 4PW 4 Lg /M eS /12 i2ro 11/3 ^iT4(wYK) 9/4 - This condi- 
tion becomes easier to satisfy as the PNS cools, allowing magnetized winds to 
dominate later stages of the KH epoch for PNSs with even relatively modest B v 
and Q. NR magnetically-driven winds, in addition to being more powerful than 
spherical, thermally-driven outflows, are efficiently hoop-stress collimated along 
the PNS rotation axis [8]. The power they deposit along the poles may produce 
asymmetry in SN ejecta distinct from the shock-launching process itself. 

Strong magnetic fields and rapid rotation can also increase the out- 
flow's power through enhanced mass-loss because E™ g oc Mq 3 . When 
the PNS's hydrostatic atmosphere is forced to co-rotate to the outflow's 
sonic radius cu s = (GMsin^LCFL]/^ 2 ) 1 ^ 3 then Mq is enhanced by a factor 
c f ~ expKvfcv/c^,,) 2 ] over M th /47i due to centrifugal ("cf") slinging [9], where 

Vfa ~ -R,,f2sin[#LCFL] ~ R V Q\J R u / ujy is the PNS rotation speed at the base 
of the open flux. Using our estimate for c s v from § 1, we see that enhanced 

mass loss becomes important for P ms < P cf , ms = ^Vo^^'KM)" 1 / 3 
(i.e., only for PNSs with considerable rotational energy E mt > 10 52 ergs). 
Fully enhanced mass loss (Mq = ^th</>cf/47r) requires u>a > cj s , which 
in turn requires that B l5 > 5 cf>15 = P^ s 4 R^ 3/4 M^ /2 _ 4 M^(uj y /uj l ) 5 / 4 ~ 
0.3P r ^ 4 Lsf e^ 3 M- 4 1/2 P^ where we have taken 

M t h from § 1. For cases with B v < B ci but P < P c f, Mq lies somewhere between 
M th /4ir and </> c fM th /47r (see [10] for numerical results). Millisecond proto-magnetars 
generally attain c f, except perhaps at early times when the PNS is quite hot. 



2.2. Relativistic Winds and Gamma-Ray Bursts 



As the PNS cools, eventually uja — > and the PNS outflow becomes relativistic 
(REL). This transition occurs after tkh for most PNSs (they become pulsars), but 
rapidly rotating proto-magnetar winds become relativistic during the KH epoch 
itself. Similar to normal pulsars, PNSs of this type lose energy at the force-free, 
"vacuum dipole" rate: E^! ~ 6 x 10 49 5f 5 P~ s 4 i^ (u;Y/^L)~ 2 ergs/s (again modulo 

corrections for excess open magnetic flux -E^ag K ^lopen ^ (^y/^l)~ 2 [8]), which 
gives a familiar spin-down timescale tj = E rot /E^^ w 300-B 15 2 P ms -Rf 4 Mi. 4 (u;Y/^L) 2 
s. On the other hand, the mass loading on a PNS's open magnetic flux is set by 
neutrino heating, a process totally different from the way that matter is extracted 
from a normal pulsar's surface. In fact, a proto-magnetar outflow's energy-to-mass 
ratio a is given by 

-arm 1 <« 

From equation ([6]) we see that because a PNS's mass- loss rate drops so precipitously 
as it cools, a oc L~ 5 / 3 e~ 10 / 3 rises rapidly with time, easily reaching ~ 10 — 1000 
during the KH epoch for typical magnetar parameters [9], [10]. Detailed evolution 
calculations indicate that E TOt is extracted roughly uniformly in log(cr) [10]. 

To conclude with a concrete example, consider a proto-magnetar with B v = 10 16 
G and P ms — 3 at t — 10 seconds after core collapse. From the cooling calculations 
of ref [3] we have L 52 (10 s) » 0.1 and eio(10 s) « 1 (see Figs. [14] and [18]) and so, 
under the conservative estimate that u Y = co> L , equation gives a ~ 500. Because a 
represents the potential Lorentz factor of the outflow (assuming efficient conversion 
of magnetic to kinetic energy), we observe that millisecond proto-magnetar birth 
provides the right mass-loading to explain gamma-ray bursts (GRBs). Further, the 
power at t = 10 s is still « 10 50 erg/s with a spin-down time tj ~ 30 s, both 

reasonable values to explain typical luminosities and durations, respectively, of 
long-duration GRBs. Lastly, because a rises so rapidly with time as the PNS cools, 
in the context of GRB internal shock models a cooling proto-magnetar outflow's 
kinetic-to-7-ray efficiency can be quite high; our calculations indicate that values 
of 10 — 50% are plausible. We conclude that magnetar birth accompanied by rapid 
rotation (but requiring less angular momentum than collapsar models) represents 
a viable long-duration GRB central engine. 
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